FTUAM-97/18 



Yang-Mills Fields on the 4-dimensional torus. 
Part I: Classical Theory [] 

A. Gonzalez-Arroyo 

Depto. Fisica Teorica C-XI, 
Universidad Autonoma de Madrid, 
Cantoblanco, Madrid 28049 , SPAIN 



Abstract We review some of the most important results obtained over 
the years on the study of Yang-Mills fields on the four dimensional torus at 
the classical level. 



1 This contains the first part of the lectures given by the author at the 1997 Advanced 
School on Non-perturbative Quantum Physics (Pehiscola) 



1 



Contents 



1 Introduction 



2 Notation and general formulas 



3 Boundary Conditions and Twist 
3 Twist Matrices 



471 Abelian twist matrices 



4.2 Twist eaters 



4.3 The algebra of twist-eaters 



5 Topology of gauge fields 



|6 from U(N) to SUjNjJZ 



N 



7 Gauge invariant Quantities 



|8 Classical Solutions! 



871 Zero-action solutions 



8.1.1 Irreducible twists 



8.1.2 Reducible twists 



8.2 Abelian solutions 



8.3 Non-abelian constant field-strength solutions 



8.4 General non-abelian minimum action solutions 



18.4.1 SU(2 



8.4.2 SU(N) 



8.4.3 SU(oo 



875 Other results on classical solutions 



|9 Fluctuations around classical solutions 



10 Concluding remarks 



1 Introduction 



Yang-Mills fields play a crucial role in our understanding of Particle Physics. 
Under a relatively simple and elegant formulation, a rich behaviour hides. 
Phenomena which are of great importance in the real world, such as Confine- 
ment and chiral symmetry breaking, are supposed to follow from the dynam- 
ics of these fields. However, these phenomena are highly non-perturbative 
and still defy ourselves to try to derive them properly from first principles. In 
trying to address such complicated issues a good deal of new ideas and strate- 
gies have been put forward by different authors over the years. Some of these 
techniques constitute a major breakthrough in our capacity to derive quan- 
titative and qualitative approaches to Quantum Field Theory (QFT). What 
in our opinion is the most important new method which arose in this context 
was the formulation of Yang-Mills theories and other QFTs on a space-time 
lattice Q. This method, not only gave rise to new non-perturbative methods 
of computation, but also provided a new illuminating way of understanding 
the phenomena taking place in Quantum Field Theory. Other major con- 
tribution followed the work of Polyakov, who argued about the relevance 
of classical configurations in computing and explaining some effects in QFT 
P, [3[. Most notably the so-called pseudo-particles, of which the best known 
example is the instanton are known to play an important role in some 
statistical mechanical systems and are expected to do the appropriate in 
QFT. 

In this review we will study a different strategy put forward by 't Hooft || 
H, [7J. He considered Yang-Mills fields living in an Euclidean space-time torus 
T 4 . In 't Hooft's original presentation, the interest of this approach lies in 
the characterization of Confinement that it provides, and the information it 
supplies on the possible phases of Yang-Mills theory. In very brief terms, we 
can say that the topology of the torus gives rise to a non-trivial topology in 
the space of Yang-Mills fields which has an appealing physical interpretation. 
The way we see it nowadays, this is just one aspect of the interest that a 
formulation of Yang-Mills fields on the torus can have. First of all, the torus 
provides a gauge invariant infrared cut-off of the theory. For that purpose, 
we could have taken any other Riemannian 4-dimensional compact manifold, 
but the torus has several advantages which we will now list: 

• The group of traslations is abelian, and hence one might still use Fourier 
decompositions. Furthermore, one can use a flat metric in all points of 
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space-time. 



• One can regard all configurations on the torus as configurations in R 4 
which are periodic. In other words, one can glue up configurations 
defined on the torus, to make up configurations on R 4 . 

• In the lattice formulation of gauge theories, numerical methods demand 
a finite number of degrees of freedom. Henceforth, one is forced to 
put the system in a finite volume. For reasons similar to the ones 
mentioned before, one usually imposes periodic boundary conditions, 
which is equivalent to having the system living on the torus. The study 
of Yang-Mills on the torus allows us to understand some of the observed 
finite size effects. One can take the continuum limit, in a way such that 
the size of the torus stays finite in physical units. 

• Topology of the torus is non-trivial (both 7Ti(T 4 ) and ifj(T 4 ) are non- 
zero). This allows the topological characterization of some properties of 
the theory, including Confinement. This is basically 't Hooft's original 
motivation. In addition, as is the case in other theories, boundary 
conditions serve to stabilize several classical configurations. 

• A curious phenomenon takes place in the large N limit. Eguchi and 
Kawai || argued, that in that limit, some quantities might have no 
finite-size effects. Hence, the torus would give results which are identi- 
cal to those of R 4 . 

• The size of the torus provides a parameter, which can be used to in- 
terpolate between the perturbative (small size) region and the non- 
perturbative (big size) one. This can provide a method for approx- 
imately calculating some non-perturbative quantities, as argued by 
Liischer ||. In addition, if the ideas of the author are correct, this 
can provide an insight into the structure of the Yang-Mills vacuum and 
the origin of Confinement [[II]] . 

To conclude, we can say that the study of Yang-Mills on the 4-D torus serves 
as an example of the behaviour of the fields on other non-simply connected 
manifolds. The latter can be relevant in other contexts. 

In these lectures, we will review the most important results obtained by 
different authors on this topic. They can be split into two parts: classical and 
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quantum theory. In this paper we will present the first part. The quantum 
part, which will follow, will make use of the results presented here. 

Our aim is to make a pedagogical and self-contained introduction to the 
subject. Accordingly, we will present all the ideas and results in the most 
simple language. We will nevertheless frequently comment upon the more 
precise mathematical statement of those results. In other cases, we will refer 
the readers to the relevant publications. 

The layout of the paper is shown in the Contents. 

2 Notation and general formulas 

Let us consider a 4-dimensional torus T 4 of size l x l x x l 2 x l 3 . The metric 
is Euclidean and we take coordinates x (x^ for /i = ... 3) on the torus as 
follows: 

< x M < Z„ . (1) 

Let us introduce the symbol \x = (0, . . . , 1^ 0, . . .) which has all but the /ith 
component equal to zero. 

We now consider Yang-Mills fields defined on T4. For simplicity we will 
restrict ourselves to the SU(N) gauge groups. We will express the vector 
potential fields in the usual matrix form: 

A,(x) = Al{x) \ a , (2) 

where A a are the generators of the SU(N) Lie algebra in the fundamental 
representation, normalized by: 

Tr(X a X b ) = 5 f 

[A a , Aft] = ifabc^c ■ (3) 

The expression for the field tensor F IMU (x) in terms of the vector potential 
A M (x) is given by: 

F^ix) = d^A u (x) - d v A^(x) - 1 [A„(ar), A„(a;)] . (4) 

Gauge transformations operate on the vector potential A /i (x) as follows: 

[q(x)]a m (x) = n(x)A fl (x)n + (x) + tn(x)d fl n + (x) . (5) 



5 



In a more precise mathematical language, Yang-Mills fields on the torus 
are connections on a principal fiber bundle whose base space is T 4 . A good 
introduction to this formulation aimed to an audience of physicists can be 
found in Ref . JTl] . In order to formulate the theory, one has first to define an 
SU (N) principal fiber bundle over T4. This can be done by giving a covering 
of the torus by open contractible sets U a . These sets can be considered 
gauge-coordinate patches. In the overlapping regions ( x 6 U a fl Up ), the 
consistency of the two descriptions implies that the two gauge-coordinates 
are related by an SU(N) gauge transformation Q a p(x), which is called a 
transition function. The open sets and the transition functions define the 
principal bundle. 

Connections define a way to parallel transport from one point of space 
to another along some path. Within one patch U a the parallel transport 
matrices are given by the well-known ordered exponential formulas in terms 
of the vector potentials A^\x) for this patch: 

U( 1 :P^Q)=T exp (-1 j A£*> {x)dx)j , (6) 

where ordering within the path goes from left to right. In the overlapping 
regions between two patches ( x G U a C\Up ), one might use instead the vector 
potentials A^\x), related to the previous ones by a gauge transformation 
involving the transition matrices: 

Af (x) = W!)]AW(i) , (7) 

which can be considered a change of coordinates for the gauge fields. The 
link with our original presentation follows by identifying our vector potential 
A M (x) with that of a single patch which covers the whole torus. 

3 Boundary Conditions and Twist 

By continuity one can extend A M (x) and F^x) to the boundary of the 
hypercube [0, Z ] x [0,/i] x [0, Z2] x [0, 13]. However, one must realize that for 
a point x located on the face x u = 0, there is a corresponding point x + z> 
which labels the same point on the torus. Henceforth, the gauge fields A^(x) 
and Ap(x + u) have to be physically equivalent. For gauge fields, equivalence 
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does not imply equal vector potentials. It is enough if they are equal modulo 
a gauge transformation: 

A M (x +u) = [Q^A^x) , (8) 

where Q u (x) are elements of the SU(N) group which depend only on the 
transverse coordinates x p with p ^ v. These matrices are known as twist 
matrices. 

If we now consider a point x located at the two-dimensional surface x v = 
x p = 0, consistency of the previous boundary conditions demands: 

A M (z> +p) = [Q u (x + p)Q p (x)]A p (0) = [Q p (x + i>)n u (x)] A M (0) . (9) 

This implies 

Q p (x + i>)Q u {x) = z pv Q u {x + p)£l p (x) , (10) 

where the constants phases which can occur due to the quadratic 

way in which Q enters in the gauge transformation formula Eq. |5|. Since 
Q p (x + u)Q u (x) and Q u (x + p)Q p (x) belong to SU(N), the phases z pv are 
elements of the center of this group (Z;y). Hence, we can express them as 

Zp U = exp{2m^} , (11) 

where n pv is an antisymmetric tensor of integers defined modulo N. This 
tensor has 6 independent coefficients, which can be expressed as usual in 
terms of two independent three vectors (k and m) as follows: 

n i = fa . (12) 

We will refer to these boundary conditions as twisted boundary conditions, 
and to n pv as the twist tensor. The elements of the tensor n M ^, being integers, 
label topologically inequivalent sectors. The particular case n^ u = is often 
referred as no-twist. In particular, choosing the twist matrices equal to the 
identity, corresponds to this case. The gauge fields are then said to be strictly 
periodic. 

Under gauge transformations, the vector potentials transform as shown 
in Eq. [5] , but the twist matrices also change as follows: 

n p (x) — ► n'^x) = n(x + p) n„(x) n + (x) . (13) 
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However, notice that the change of twist matrices does not produce any 
change of the twist tensor n^ v . 

In addition to gauge transformations, there is another symmetry group 
present. It is the group of transformations: 

%{x) — > a^x) = z» %{x) , (14) 

where z M are elements of the center. This group is isomorphic to Z A N and 
plays an important role in what follows. 

The previous introduction of twisted boundary conditions follows very 
closely the one made by 't Hooft in his first papers on the subject |5|. 
Using the language of connections on fiber bundles, the twist matrices are 
simply compositions of different transition functions in the limit when one of 
the patches covers the whole torus | I2| , |T3|| . In this respect, the consistency 



condition Eq. ^ is related to the so-called co-cycle condition, that transition 
matrices have to satisfy: 

Q al = Q a( gf2g 7 . (15) 

The occurence of the phases z^ u and of the twist tensor, has to do with 
the fact that gauge transformations of the vector potentials are insensitive 
to the center of the group Z^. Henceforth, one could consider the gauge 
group (structure group) as being SU(N)/Zn. A non-zero twist tensor can be 
seen as an obstruction to go from an SU(N) /Z^ principal fiber bundle to an 
SU(N) bundle. This obstruction can be related to the second homology class 
of the base manifold with coefficients in the center of the group H 2 {T^ Z^). 



For a proof of these facts see for example Ref. [14 . 



4 Twist Matrices 

The first problem which we will address is the form and existence of solutions 



to the twist matrix equations Eq. [10]. We will consider two classes of solutions 
which will be relevant in what follows. The first class is given by abelian 
space-dependent matrices. The second class is given by constant non-abelian 
matrices: twist eaters. 

4.1 Abelian twist matrices 

Here we will show that there is always a solution to the problem of finding 
twist matrices satisfying Eq. [1(1. This solution is given in terms of commuting 



8 



matrices, and is hence referred as abelian. Let us set: 

Qp(x) = exp {% uj^(x)} , (16) 

with [uj^{x) , lu v (x) } = 0. If we diagonalize the matrices u^, the twist condi- 
tion Eq. [10] amounts to the following condition on its eigenvalues a;": 

077-77 

A u u;(x) - A^x) = + 2tt<£, , (17) 

where Q M „ =diag(q^,) is an integer diagonal matrix satisfying TV(Q^) = 
—n^, and the symbol A M is defined as: 

A^{x) = 4>{x + p) - <j>{x) . (18) 

A particular solution of the previous equation is given by: 



with T = diag(l, . . . , 1, 1 — N). This solution corresponds to the choice 
Q^iv = for a 7^ and q^ u = —n^ u . For arbitrary value of the integers q^ v 
one gets: 

= , (20) 

with 

TV = n^l + NQ,j, v . (21) 

A general solution can be obtained by adding to the previous solutions a term 
A^<ft(x), which is a general solution of the homogenous equation. This general 
solution can be easily seen to be a gauge transformation of the solution 



expressed in Eqs. ETJ- 21 



4.2 Twist eaters 

A class of solutions which plays an important role in the following are the 
so-called twist eaters or twist eating solutions. These are constant matrices 
— satisfying: 

r^ = ex P {^^}rjv (22) 
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The previous equation defines a group, whose generators are T M and the 
phase exp {2m /N}, and satisfying the relations Eq. ^ involving {p^ u = 
We will call this group the twist group Q . Notice that the Clifford Algebra 
is a particular case of this class of groups. Our solutions can be looked at 
as unitary representations of this group. Hence, we can use the standard 
methods and ideas on representations. We will actually be more general that 
required and work in arbitrary space-time dimension d (jU = . . . (d — 1)). 

Given one solution T M one can construct another one by means of: 

r M — > r; = fiiy^ + (23) 
r M — r; = w, (24) 

with Q G SU(N) and G Zat. Furthermore, given two solutions with 

n W (3) 

W = Ivf = Iv 1 ' one can cons t ruc f the direct sum, which is an Ni + N 2 
dimensional representation. The new representation is reducible. One needs 
only to find the irreducible representations and the rest can be formed by 
direct sum. 

The choice of generators is not unique, and thus different sets of (p^ 
correspond to the same group. To investigate this point let us introduce the 
following notation for the product of the generators: 

r(n) = rs»i? ...r^\ (25) 

where are integers. The matrices T(n) satisfy: 

r(n)r(ra) = T(m + n) exp{ — ^ m v n^} . (26) 

The matrices of the form zT(n) form the group Q. One might choose other 
generators r' = r(s^) provided the matrix s£ is invertible. Notice that if 
we use r' instead of T M , the twist tensor changes to: 

— = EaW s^n pa . (27) 

For irreducible representations, by Schur lemma, all the matrices which con- 
mute with the generators (elements of the center) must be multiples of the 
identity. In particular, any matrix of the group taken to the Nth power 
belongs to the center. This allows the numbers s„ to be considered inte- 
gers modiV. This enlarges the group of transformations of generators from 
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SL(d,Z) to GL(d,Z N ). One can use these transformations to bring ip^ u to 
a canonical form: 



if^ u (canonical) 



A 

-A 





V u u / 



o 





(28) 



where A = diag(^-, ...,—). The integers Pi > 1 are such that pi divides pj for 
i > j. The value of r < [|] will be referred to as the rank of the group Q. An 
important subgroup is the center C of Q. It is generated by (r' i _ 1 ) p % (T' r+i _ 1 ) Pi 
for i = l...r and T' 2r . . . r' cZ _ 1 . For irreducible representations, all these 
matrices have to be a multiple of the identity r M I, and the constants 
characterize the representation in question. 

Now let us use the standard strategy of considering the maximal abelian 
subgroup A. It is generated by , i — 1 . . . r (factoring out the center). 
Let | A > denote a vector which is a simultaneous eigenstate of these matrices 
with eigenvalues Aj_i (which can be written as an r dimensional weight vector 
A). The matrices T^, +fc _ 1 now act on this state and transform it into another 
eigenstate with eigenvalue (t k X)i = Aj_i exp{^5 ifc }. In this way, given one 

state |A > one generates P1P2 ■■■p r different states |A > ( Notice that 
Aj-i 7^ since it is a representation of a group, and the matrices have to be 
invertible). Now we can normalize these eigenstates so that: 

r;. +fc -i|A>=x fc |t fc A> . (29) 

The definition is consistent because the T' r+k _ 1 with different k conmute. We 
have to insert Xk because \ P k = T r+k-i is a Casimir, which is not necessarily 
equal to 1. 

Thus, we have explicitly constructed the irreducible representations of 
the twist group. The dimension of the representation is: 

N = pi x p 2 x ■ ■ -p r . (30) 

Henceforth, supposing that there is a solution in N x N matrices, then the 
representation is irreducible if and only if iV = iVo- One can construct 
reducible solutions in all dimensions which are multiples of Nq. 

The irreducible representation is unique modulo similarity transforma- 



tions (Eq. |23|) and multiplication of the matrices by a constant (Eq. 24) 
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Representations that are related by a similarity transformation are called 
equivalent. As we can see, equivalent representations are labelled by the 
Casimirs: Tj_i = Afl 1; r r+i _i = xT (f° r i — and r M (for /x = 

2r, . . . , d — 1). For unitary representations these Casimirs are phases. We 
are mainly interested in representations in terms of SU (N) matrices. In that 
case, the phases have a discrete set of values (e Z N ) and the number of in- 
equivalent irreducible representations is discrete. We will now calculate how 
many of them there are, since this will turn out to be relevant later on. We 
get: 

#ofirr. reps. = (^-J . . . N*-* = K~ 2 ■ (31) 

This follows from counting the possible values of the Casimirs consistent with 
the requirement that the matrices belong to SU(N). Having studied the 
problem in general, let us now take a brief look at the 2, 3 and 4-dimensional 
cases, which are the interesting ones for us. 

In 2 dimensions, we have 2 matrices T and r x and a single twist tensor 
element uq\ = (foiN. If q — gcd(n i, N) then there 3 a couple of integers n\ 
and ri2 such that: 

n±N + n 2 n 01 = q . (32) 

Now the canonical generators are T' = Tq 2 , T[ = F±, and p\ = y. The 
representation is irreducible if pi = N (=>■ q — 1). Our canonical basis 
amounts to r' = and T[ = where: 

(Q N ) kj = z6 jk exp{^} , (33) 

where indices are defined modulo N, and z is chosen to ensure determinant 
equal to unity (z=l for odd N and exp{^} for even N). Thus, the solution to 
our problem is T = Q^ 01 and r x = P^, which is unique modulo similarity 
transformations (if N and n i are coprime). 

In three dimensions we have n^- = ei^m k . Then the 3 canonical matrices 
are T[ = r(s^ 1 ^), T' 2 = r(s^ 2 ))and r' 3 = T(m/mo), where m = gcd(m, N). 
The vector can be chosen among those that satisfy: gcd(s^ 1 ^ ) x m, N) = 
m . The other vector is such that (s^ 1 ^ x m) ■ ^ = m mod N. One can 
easily see that p\ = N/m and the representation is irreducible for m = 1. 
In this case the generator of the center of the group must be a multiple of the 
identity: T' 3 = z 3 I. The inequivalent irreducible representations are labelled 
by z 3 . 
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Finally, we arrive at the 4- dimensional case. We know that there are 2 
subcases: rank 1 and rank 2 (rank is trivial). The relevant quantities are 
Pi and P2 (for rank 2). In order to obtain these numbers it is not necessary 
to bring n^ v (</V/) to the canonical form. They can be determined in terms 
of invariants of the transformations. We see that q\ = gcd(n^, iV) is an 
invariant. By looking at the canonical form, one sees that: 

N . . 

p ' = Sdf^i • (34) 

If we perform SL{4, Z) transformations to n^ U) there is another invariant: 



fJ,U J 



-n^n^ = k-m = Pf{n^) , (35) 
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where n Mi , = (1/2) e M ^ p(J n pa and Pfin^) stands for the Pfaffian of the an- 
tisymmetric matrix n^ v . However, n{n^ v ) is non-invariant under the wider 
group of transformations GL(4, Zjy), or under changes of by integer mul- 
tiples of N. By first bringing n pu to a canonical form with SL(4, Z) transfor- 
mations and later applying the wider class of transformations, it is not hard 
to see that the other invariant P2 can be obtained as follows: 

N N 

P2 = 7 s = • (36) 

We remind you that the necessary and sufficient condition for the existence 
of solutions is that iVo = pip% divides N . The representation is irreducible 
for Nq = N. The case of rank 1 corresponds to p% = 1 4^ qi = N. One can 
summarise our results in the following statements: 

• The necessary and sufficient condition for the existence of solutions is 
that gcd(«(n AtI/ ), N) = N. This case is called Orthogonal twist, since 
it corresponds to k ■ rh = mod N. The opposite case is refered to as 
Non-Orthogonal twist. 

• The representation is irreducible, and hence unique modulo the trans- 
formations |2|]-|23], provided gcd(iV, n^, «(n^)/iV) = 1. It must be one 



of the N 2 inequivalent irreducible representations. 
The representation has rank 2 provided gcd(^y^, N) ^ N 
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Let us conclude this section by pointing to the relevant references where 
the study of twist eating solutions was addressed. The first solutions of this 
type were found in Ref ]I5[ and jn|. 't Hooft addressed and solved the 





This also covers the 2 
IfJ. Rank 2 



problem for the case of rank=l in 4 dimensions 
and 3 dimensional cases, the latter studied previously in Ref. 
solutions, whose existence was established in Ref. [Q, appear and become 
relevant in the Twisted- Eguchi-Kawai model |18|, [19| . This triggered renewed 
interest in finding a complete solution to the problem. The final solution up 
to 4 dimensions was obtained in Refs. |20], |2T|. Finally, the d-dimensional 
case was studied in Refs. H22|, |23 



4.3 The algebra of twist-eaters 

In this subsection we would like to consider the case of irreducible twists. 
Let us examine here what would be the equivalent of the Clifford algebra 
for our case. We consider the T(n) matrices defined in Eq. We will see 
in what follows that there are as many linearly independent such matrices 
as there are elements in the quotient group Q/C. This group is abelian and 
isomorphic to Z P1 x Z P1 x . . . x Z Pr x Z Pr . As a consequence of our study of 
the last subsection we know that the order of this group is iV 2 . 

Since the representation is irreducible the matrices corresponding to an 
element of the center group C are multiples of the identity. By virtue of the 
multiplication formula Eq. the matrices corresponding to two elements 
belonging to the same class in Q/C, differ by multiplication by a phase. 
Hence, at most there are as many linearly independent matrices as elements 
in Q/C. Now let us prove the following result: 

Theorem 1 If for each class a 6 Q/C we choose a representative: 

A( a ) = e^ (a) r(n(a)) , (37) 
the corresponding matrices are linearly independent. 

The proof is very similar to the corresponding one for the Clifford algebra. 
One makes use of the following result: 

Lemma 1 Let A and B be two invertible N x N matrices such that A B = 
zBA with Zy^l, then z 6 Z N and Tr(A) = Tr(B) = 0. 
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This lemma can be easily proven. By diagonalizing A, for example, one 
reaches the condition (An — z Ajj) Bij = 0. Now since A is invertible An ^ 0, 
and hence Bn = 0. This shows that indeed the trace of B is zero. Now, since 
B is invertible, for any i there must exist i' such that Bn> ^ 0. Hence, for 
this pair we must have An — z Aw = 0. If we start with i' and repeat the 
argument, there must exist another i" satisfying the same relation. Since the 
number of different eigenvalues of A is finite, this implies that there is an 
integer p such that z~ p = 1. Hence starting from one eigenvalue one gen- 
erates a sequence of p different ones. Now suppose A has one eigenvalue r 
with degeneracy s. Then as we have seen z~ x r must also be an eigenvalue. 
In addition, its degeneracy must be s as well. This is so because B induces 
an homomorphism between the spaces of both eigenvalues, and due to the 
invertibility of B, it must be an isomorphism. In this way, from each eigen- 
value one generates ps vectors. Repeating the operation for an eigenvalue 
not contained in this set, we conclude N = p Sj), and hence p divides N. 
By diagonalising B instead, one can show that the trace of A is also zero. 
This completes the proof of the lemma. In fact, our lemma is very much 
related to the explicit construction of irreducible representations of Q done 
in the previous subsection. 

As a consequence, the trace of all T(n) matrices is zero, except for those 
which are elements of the center group. Now let us proceed to prove the 
linear independence, asserted in the Theorem. We will work by reductio ad 
adsurdum. Consider that there exist a non-trivial linear combination of the 
A (a) matrices which is zero: 

c{a) A(a) = . (38) 

aeg/c 

Then, if we multiply both sides of the equation by r _1 (ri(a)) and take the 
trace, we obtain N e l ^ a ^c a = 0. Since this is true for all a, all the coefficients 
must be zero, thus contradicting the non-triviality assumption. QED 

Our previous theorem impies that the A (a) matrices define a basis of the 
complex vector space of all N x N matrices. The product of these matrices 
satisfies: 

A(a)A(fo) = e ld{a > b) \{a + b) , (39) 
where, by virtue of Eq. ^6], d(a, b) is given by: 

d(a, b) = 0(a) + <p{b) — <p(a + b) + n^a) n u (b) mod 2n . (40) 
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The phases d(a,b) depend upon the choice of representative (4>(a), n(a)). 
However, the antisymmetric combination: 



< a, b >= d(a, b) — d(b, a) = n M „ n^(a) n u (b) mod 2n (41) 



fJbV 



does not depend on the choice. 

The real subspace of hermitian N x N matrices is isomorphic to the Lie 
Algebra of U(N). This space can be obtained by linear combinations of the 
A (a) matrices with coefficients satisfying certain relations. Using the freedom 
to make specific choices of the representatives (0(a), n(a)), one can simplify 
the form of the relations. For example, by choosing the matrices, such that 
A(— a) = A + (a), the condition of hermiticity reads c(— a) = c*(a). Hence, one 
can obtain the form of the structure constants of U(N) in the basis of A (a) 
matrices, as follows: 

fabc = -t5ca +b (e ld{a > b) -e^) . (42) 

Other choices of representatives can impose certain properties on the A (a) 
matrices. 



The algebra of twist-eaters was introduced in Refs. [p4|, |19|j in 2 and 
4 dimensions respectively. In those references, it was used as a basis of 
the Lie Algebra of SU(N). This will turn out to be useful later (See also 



Refs. 25 



5 Topology of gauge fields 

In addition to twist, we must consider the ordinary topological charge (in- 
stanton number): 

Q = TT^ I T < F ^( X ) *V(z)) d A x . (43) 
itOT Jt 4 

Using the fact that Tr(F AiJy (x) F^(x)) is a pure divergence (d^IC 1 ), and in- 
tegrating once, we get: 

<2 = £/^ A ,*'*. («) 
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where A p is defined in Eq. [T3|, da p is the integral over the 3 dimensional face 
= 0, and is given by: 

K» = 2 e^ vpa Tr(A u d p A a - ~z A v A p A a ) . (45) 

The main ingredient that we need now, is how does change under a gauge 
transformation. We have: 

K"([n]A) = K»(A) - 2 -% Tr(A v A p A a ) - 2 e pupa d p Tr(A u A a ) , (46) 
with A u = iVL + d v Vt. In our case A v (fl) = [Qp]A u (0), and hence we get: 

Q= -^e^jda^A^A^A^) (47) 
— s~5 e nvpcr J da^p Tr(A p (AJ* A a )) , 

with da pp the integral over the 2-dimensional surface x p = x p = 0, and 
A^ ] =iVL+d u VL p . Now using D P A^ = D P A { U P \ with D p defined by 

DpA^ = n+(A^) + AM)n p - AS* , 

and manipulating the previous equation, we arrive at: 

Q= -^e p ^jda p Tr(n+(d u n p )n+(d p n p )n;(dM,)) (48) 



/ da pp Tr((A p 4 M) + AY>) 



8tt 2 c /^P°" 



which is our final formula. As expected, only the twist matrices Q p enter 
the formula. Hence, given the twist matrices, one can compute the topo- 
logical charge. Since for every twist we have found a particular set of twist 
matrices — the abelian ones — , we might compute the value of Q for this class 



of solutions. By taking Eqs. |iq , [20| , pi] and substituting it into Eq. [48], we 
obtain: 

Q = + \Tr{QM = + . (49) 

The second term is an integer. For example, had we taken the particular 
solution Eq. [H| we would have gotten Q = K{n pv )^-^-. Actually, with a 
suitable choice of e« = Qoi and f3 i = jQjk^ijk, one can make the second term 
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of Eq. take any possible integer value n. 
job is (i = 1, 2): 



A possible choice that does the 
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The previous formula makes sense for N > 3. For SU(2) and no-twist 
(k = m = mod 2), there is an exception to the rule: Only even values of 
the topological charge are attainable in this case. 

In summary, we have shown that in the presence of twist, the topological 
charge is given by: 

Q = -^%^ + n, neZ. (51) 



N 

The previous formula was conjectured by 't Hooft |5|, |J, and actually proven 
along similar lines to our presentation by Pierre van Baal (12| . A consequence 
one can derive from Eq. |5l| is that for non-orthogonal twists (^(n^) ^ mod 



N) Q is not an integer. For orthogonal twists, one has twist-eating solutions, 
whose topological charge vanishes. 

Now let us comment about the appropriate mathematical explanation 
of our results. The topological charge for an SU(N) group, is just minus 
the second chern number, which is an integer. In view of this, the previous 
result is perplexing. However, as mentioned previously, we have to look at 
our fields and bundles as being defined in SU(N)/Zn. Hence, we should 
rather compute the previous expressions in the adjoint representation, which 
is a faithful representation of SU(N)/Zn. To do so, we have to replace in 
Eq. H A a by X^^l The result gets multiplied by = 2N, where c a ,f 
and d,A,F are the quadratic Casimir and the dimension of the adjoint and 
fundamental representations. We get: 



Q 



(adjoint) 



2k(?v) + 2Nn . 



(52) 
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which is an integer. The previous quantity is more appropriately referred to 
as the first Pontryagin index pi. Indeed, mathematicians have proven that 
px and twist provide a complete classification of the bundles jl4 . 



6 from U(N) to SU(N)/Z 



N 



In this section, we will look at twist from a different point of view, by con- 
structing a twisted SU (N) gauge field as a projection of a U (N) gauge field 
one[fZ7]]. Our starting point, hence, is a U(N) gauge field on the torus: 
A' (x) = A"^(x) X a ■ The difference with respect to the SU(N) field case, 
is the presence of an additional generator Ao = This field satisfies the 

boundary condition: 

A'^x + P) = [tf v (x)] Afr) , (53) 

where now the matrix £l' v {x) belongs to U(N). The compatibility condition 
satisfied by these matrices is: 

n , p (x + p)n' u (x) = n' u (x + p)n' p (x) . (54) 

It is now possible to decompose the gauge field into an SU(N) part A M (x) 
and a £7(1) part, by considering the corresponding components of A'Jx). 



The boundary condition Eq. [53| somewhat mixes the two parts. However, we 



might factor the matrices £l' v {x) into an SU(N) and a £7(1) part as follows: 

Sl' u (x) = exp{zu u (x)} n v (x) , (55) 

where £l v (x) G SU(N) The decomposition is unique modulo multiplication 
of Q u (x) by an element of the center z v accompanied by the corresponding 
shift in uj v {x). However, by continuity z v is space independent, and does not 
influence the following. Now, for the 2 different parts we have the following 
boundary conditions: 

A M (x + v) = [Q u (x)]A^x) (56) 
A>l(x + P) = A'l(x) + V2Nd fl u l/ (x) . (57) 

The compatibility conditions read: 

Q p (x + P)Q v (x) = exp{2vr^} Q u {x + p)Q p {x) (58) 

A p lu u (x) - A u to p (x) = ~^ n pv m od 2 7r . (59) 
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Henceforth, in general, the resulting SU(N) field has a non-trivial twist 
tensor n pu . We can compute the field-strength tensor and decompose it into 
the SU(N) and U(l) part: 

F^(x) = F, u (x) + F' II u (x) \ (60) 

where F°,(x) = d^A'^x) - d v A'l(x). 

Now, we could use these results to relate the topological properties of the 
U(N) field to those of SU(N) gauge fields with twist. It is well-known from 
the theory of characteristic classes [II| that the topological properties of these 



fields (actually of the bundles or transition matrices) are characterised by the 
Chern classes. In our case, we have the following topological invariants: 

c^u = J dx^ A dx u Tr(F^) (61) 

/ Tr(F;,(x)F^(x))rf 4 x . (62) 
JTa 



16vr 

The first set of invariants c pu is the integral of the first Chern class over 
the non-contractible 2-dimensional surfaces of the torus. These numbers are 
integers and hence do not depend on the actual surface chosen but only on 
the directions /i and v. The other quantity k is again an integer and is 
obtained from the second Chern class. Now, we can evaluate these invariants 
in terms of our decomposed fields: 

N 

<V = — (A^u^x) - Afj,uj v (x)) (63) 
k = Q + jz7 c^Cfu, , (64) 
where we have used the boundary conditions and compatibility equations. 



Notice, that using Eq. |59] one arrives at: 

*V = <V mocl N > ( 65 ) 

which relates the twist tensor with the first Chern class of the original U (N) 
field. Given this, Eq. |63| reproduces the expression of the topological charge 
Q of an SU (N) field with twist that we found before (Eq. [51]) . 
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7 Gauge invariant Quantities 



The most important quantities of the theory are the Wilson loops. They are 
given by: 

W(j) = Tr(U(rf : P -»• P)) , (66) 

where {7(7 : P —>■ P) is the parallel transporter matrix along the closed loop 
7. Closed paths on the torus can be classified into the different homotopy 
classes, which form the group 717 (T 4 ) = Z%. An element w = (w^) of this 
group contains all paths which wind times around the /z-th direction of the 
torus. The corresponding Wilson loops are called Polyakov loops or Polyakov 
lines. The simplest of these, are straight line Polyakov loops along the /i-th 
direction: 

P^x) = Tr(U^x)) =Tr(Texp{-i [ " A^x) dx^} Q^x, x M = 0)) . (67) 

Jo 

The ocurrence of the twist matrix f2 M in the previous formula, is essential 
to guarantee gauge invariance. In general, one should insert a twist matrix 
every time that the loop 7 traverses the boundary of the patch. 

Polyakov lines are non-invariant under the discrete group of trans- 
formations (Eq. |l^). If k = (kfj) is an element of this group, the Polyakov 
loops transform ^(7) — > exp {2m k ■ w^/N} ^(7). It is also important to 
notice that Polyakov lines are non-periodic. For example, from its definition 
Eq. |67| one can easily deduce that: 

P M (x + 0) = z^P^x) . (68) 

8 Classical Solutions 

This section will try to summarize what is known about solutions of the eu- 
clidean classical equations of motion of Yang-Mills fields on the torus. Clas- 
sical solutions are extrema of the Yang-Mills (euclidean) action functional: 

S = \ ! d 4 xTr{Y^ v {x) F^(x)) . (69) 

They satisfy D lJ F^ v {x) = 0. Some of the classical solutions are those which 
minimize the action within each twist and topological charge sector, and are 
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referred to as minimum action solutions. Using the Schwartz inequality one 
can derive the following bound for S: 



The bound is saturated by self-dual and anti-self-dual solutions (F Mi ,(x) = 
±F M j,(x)), provided they exist. For non-orthogonal twists, the topological 
charge is never an integer, and henceforth the action can never be zero. 

It is clear that given one classical solution, one can obtain others for the 
same twist, topological charge and action, by making a gauge transformation. 
They are said to be gauge-equivalent. Our purpose would be to find all 
gauge-inequivalent classical solutions. In the next sections we will review the 
different cases which have been studied in the literature. 

8.1 Zero-action solutions 

As mentioned previously, zero action solutions (S = 0) can only occur for 
orthogonal twists, which is the class of twists for which k ■ m = nin^) = 
mod N. If we take a zero vector potential A At (x) = 0, then this is only 
compatible with the boundary conditions Eq. ^ provided the twist matrices 
are constant. Indeed, this is the class of twist matrices which was studied 
previously under the name of twist eaters. It is now clear what is the origin 
of the name. 

Having found zero-action solutions for all orthogonal twists, we turn now 
to the problem of determining how many gauge-inequivalent sets are there. 
Within the set of A^(x) = configurations, one is still free to make global 
gauge transformations, which will not change the value of the vector poten- 
tial. However, these operations transform the twist matrices by a similarity 
transformation. Henceforth, the set of inequivalent twist-eaters labels the 
gauge-inequivalent configurations of this type. 

Alternatively, one might take fixed twist matrices fl^x) = T^, and con- 
sider different values of the vector potential. Since S = =>- F^x) = 0, 
there must exist Q(x) such that: 



However, the vector potentials must satisfy the boundary condition Eq. |8|, 
which in our case demands: 



S > 8vr 2 |Q| . 



(70) 




(71) 



Q(x + P)d^VL + (x + 0) = Y y n(x)d fl Q + (x) T 



+ 



(72) 
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This can be shown to imply: 

fyf v {x) = , (73) 
where f „ = Q + (x + v)T v Vt(x) . (74) 

Now the question is whether this is simply a gauge transformation of the 
previous ( A At (a;) = 0) solution. If we perform a gauge transformation with 
Q + (x), we indeed get A M (x) — * ; — ► T M . The new twist matrices 
are constant as a consequence of Eq. [73|, and we are back to the previous 
situation. In summary, we have proven that: 

• The class of gauge inequivalent zero-action solutions coincides with the 
class of inequivalent twist-eaters. 

Let us now characterize the elements of this class in terms of gauge invariant 
quantities. If we evaluate all Wilson loops for the configuration with zero 
vector potential and twist matrices equal to T M , we get 

Wiri) = Tr(r£ • • • r* ) = e«*» Tr(r(w( 7 ))) , (75) 

where s is the number of times the path 7 crosses the edge of the hypercube 
FLIP) ^t] an d a (7) a path-dependent phase. Every time the curve 7 crosses 
the edge of the hypercube along the /ith direction one gets a factor of 
in the previous expression, with e = ±1 depending on the sense of crossing. 
The final trace is only non-zero when r(u;(7)) is an element of the center C 
of the twist group. Henceforth, the Casimirs of the twist group Q correspond 
to non-zero Polyakov lines. 

We have set the mathematical problem properly: the space of gauge in- 
equivalent zero-action solutions (or euclidean vacua), which is known as the 
euclidean vacuum valley, is the space of equivalence classes of N-dimensional 
representations of the twist group (In the mathematical literature, zero- 
action solutions are referred to as flat connections). In what follows, we 
will study the structure of this space for the different twists. The study 
naturally splits into the reducible and irreducible twist cases. 

8.1.1 Irreducible twists 

Under this name we include the values of n^ v and N for which we have 
an irreducible representation of the twist group. We recall that this occurs 
when pi ■ p 2 = N for rank 2, or pi = N for rank 1. In either case we proved 
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previously that the number of inequivalent solutions is discrete. Indeed, we 
showed that the number of such solutions is iV 2 . There is always a Polyakov 
line which takes a different value in any two of these vacua. As a particular 
example, consider the twist m = (1,1,1); A; = (0,0,0). This has rank 1 
and pi = N. The two independent Polyakov lines are P = Tr(U )/N and 
P (111) = Tr(U 1 U 2 U 3 )/N. Each one can take N different values ( the N th 
roots of unity). 

Another example of rank 2 is m = (n, 0,0); fc = (-,0,0) (n divides N). 

In this case the independent Polyakov lines are: Tr([/^ n )/Af, Tr(U^ n )/N, 
Tt{Uq)/N and Tr{U™)/N . The first 2 can take n values and the last two 
N/n values. 

8.1.2 Reducible twists 

This case is much more complicated than the previous one. We have N = 
sN , where s is the number of irreducible components. One can bring the 
matrices to block diagonal form. In each Nq x Nq box one has an irreducible 
representation of the twist group. Since unitary irreducible representations 
are unique modulo multiplication by a phase, we might write all the matrices 
in the form: 

r M = s, ® rf»> , (76) 

where T^ ^ are the matrices of one irreducible representation, and are 
diagonal SU(s) matrices. 

The case that we will study first is the purely periodic one {n^ v = 0). This 
is a particular case of the general reducible case having s = N and r M = S^. It 
might seem that the eigenvalues (elements) of label the different non-gauge 
equivalent solutions. In that case, the vacuum valley would be isomorphic 
with Sf N , where d=4 is the space-time dimension. However, this is not so, 
because there are different points of this manifold which are gauge equivalent. 
This occurs because there are similarity transformations which bring the 

back to diagonal form, namely those which exchange the order of the 
eigenvalues simultaneously in all d matrices. An additional complication 
comes from the requirement that the matrices have determinant equal to 1. 
This conditions fixes the value of the last eigenvalue in terms of the previous 
ones. Hence, we have the manifold Sf^ N *\ Now we may introduce the 
group of transformations T in this manifold corresponding to an exchange 
of the N eigenvalues. This group is isomorphic to the permutation group of 
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N elements. Hence our vacuum valley is: 

V = Sf^/T . (77) 

The disgusting feature of this space, is that it is an orbifold and not a man- 
ifold, corresponding to the fact that there are fixed points of the group of 
transformations T . They occur when there are degenerate eigenvalues for all 
S M at the same time. 

The presence of an infinite number of gauge-inequivalent vacua for pe- 



riodic gauge fields on the torus was pointed out first in Ref. |pq| . It was 
realized that there are special points in this space corresponding to the sin- 
gular points. The name toron was used to refer to a gauge-equivalence class 
of vacua. In that reference, instead of setting A M (x) = and labelling the 
torons by the twist matrices, the latter were fixed to unity, and each toron 
was labelled by the value of A M (x). The transformations which leave the 
unit twist matrices invariant are the periodic gauge transformations. It is 
possible to choose one representative of each class of gauge-equivalent vector 
potentials, within the set of constant diagonal vector potentials. This is an 
alternative and equivalent description of the vacuum valley. 

If we go to the generic irreducible case, our analysis and description re- 
mains valid after the substition N — > s. Henceforth, the euclidean vacuum 
valley V is given by: 

V ee Sf'-V/T, , (78) 
where T s is isomorphic to the permutation group of s elements. 

8.2 Abelian solutions 

The first class of solutions which were known for orthogonal and non-orthogonal 



twists are abelian. Consider first the abelian twist matrices given in Eqs. [16], ^0], |21 
If F^vix) and A^(x) are chosen to conmute with these matrices, the vector 
potential should satisfy the following boundary condition: 

AvA.ix) = d»u u = ^ . (79) 

A solution of this equation is given by: 

A,(x) = -^ = -^7^V (80) 
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This is indeed a classical solution (D fl F^ u (x) = d^F^ u (x) = 0). We can 
compute for this configuration the value of the action and the topological 
charge: 

Q = ^Tr(T flu f flu ) = ^Tr(SB) , (83) 
where T 0i = £ h = e ijk B k and 

H = —ftfT- ■ ( 84 ) 

As expected, the value of the topological charge is the one which follows 
from the twist matrices. Therefore, there are solutions for all values of the 
twist and the topological charge (excepting the odd integers for SU(2) with 
twist k = rh = mod 2). The corresponding value of the action follows 
from Eq. Here we will look in bigger detail to the self-dual (or anti- 
self-dual) solutions, which when existing have the minimum possible action 
within each twist and topological charge sector. The self-duality condition 
implies = Bi. As a consequence, the g« have to be rational numbers: 

q% = ^7 , with pi,p { e Z + , (85) 
Pi 

and pi and p\ coprime. Henceforth, we might express Ei = p\ Kj and Bi = 
Pi Kj, where Kj are traceless, integer, diagonal matrices. One gets: 

3 1 

^ = E ^PiPjTr(K?). (86) 

i=l 

We still have to put in the condition Eq. [21j. Since there exist integers aj and 
bi such that ctiPi + bip\ = 1 (i fixed), we may solve for Kj: 

Ki = Si + NSi, (87) 

where Sj are integer, diagonal matrices satisfying Tr(Sj) = — s», and < 
S{ < N is such that hi = p\ Sj mod N and rrii = pi Sj mod N. The topological 
charge becomes: 

Q = Y,PiP'i f-| + Tr (S?)) >0. (88) 
i=i ^ ' 
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The first term is up to an integer equal to — K{n ^ v> . For fixed size and twist 
the minimum value of the topological charge is Q = Y^iPiPi s «(-^ — Si)/N. If 
we are given the twist (0 < ki, mi < N), it is always possible to find a given 
size of the torus for which there is an abelian self-dual solution. 
We might summarize our results in the following statements: 

• There are abelian solutions for all sizes, twists and topological charges, 
except for SU(2), no- twist and odd Q. 

• Self-dual ( Q ^ 0) solutions occur for all twists and for some torus sizes, 
such that the ratios of the lengths in each direction are square roots 
of rational numbers: U/Iq = ^/qi^Qs/ Qi- The minimum value of the 
topological charge for this configuration is Q = £\ pip^s^N — Si)/N, 
where pi and p\ are defined in Eq. |85] and Si = Oj rrii + hi ki (with 
aiPi + b i p' i = 1). 

• The only case in which these abelian solutions are minimum action 
solutions within each twist sector, is for SU(2) and k — rh — (1,0,0) 
(pi = p'i = 1). By a misfortunate coincidence, this configuration was 
called toron by the authors of Ref. |2£| — the same name used before 
for a different object — . 

The anti-self-dual solutions have negative values of the topological charge 
and can be obtained from the previous ones by time reversal or parity 



8.3 Non-abelian constant field-strength solutions 



't Hooft [|T7| discovered that a wider class of constant conmuting F fiu (x) 
solutions are compatible with the twisted boundary conditions. To find them, 
let us work the other way round and start with the expression of the vector 
potential and of the field strength tensor Eqs. |80|-[8l| We will try to find out 
what are the possible twist matrices fi M (x) in this case. They must conmute 
with Tp V . From the boundary conditions for A M (a;) (§) one gets: 



itt v (x)dSl+{x) = A v AJx) 



vr T, 



N I 



(89) 



The general solution of the previous equation is given by: 

Q^x) = exp{z — j- T^} , 



(90) 
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where T M are constant SU(N) matrices which conmute with T^ u . These 
matrices satisfy: 

r M r,r+r+ = exp{-2^Q;j . (91) 

The form of Q' follows from our previous study of twist eaters. If there are 
Ni coinciding eigenvalues in T^, we might set T M within this subspace to 
be one of the twist eating solutions. Hence, in this subspace, one has that 
the eigenvalues of Q' are given by q'$ /Ni, with qffi an orthogonal twist 
tensor in SU(Ni). Now imposing the twist condition Eq. |10| we get: 

T AB/ = n AB/ I + JVQ^ . (92) 

The difference with respect to the purely abelian case lies in the expression of 
(versus Eq. [H]) , since Q' is not an integer matrix. With this difference, 
expressions jS2] - |8~3| are valid in this case as well. If we focus on self-dual 
solutions, the main difference with the previous case is the form of Kj, in 
which Eq. |7] is replaced by: 

(a) 

Ki = 8i + N^-6 a . (93) 

What we have done, is to separate the N eigenvalues into sets of N a coinciding 
ones, so that ^ a N a = N. The matrix 5 a is the projector onto the subspace 
a, and are integers such that = — Sj. The formula for the 

topological charge is now: 

1=1 \ a / 

The abelian formula Eq. RSI is a particular case of this one for N a = 1 Va = 
1 . . . N. The condition that the twist in the a subspace be orthogonal, 
amounts simply to the fact that the contribution of this space to the previous 

formula be an integer pip[ ^ G Z). 

The particular case studied explicitly by 't Hooft is that of 2 subspaces 
iV = iVi + iVV We can work out explicitly the topological charge for this 
case, and we get: 

Q = Ypp> {slN - SiNl)2 (95) 
^ 2^P*P* NlN2N W 

i=l 
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In this shown by 't Hooft ]T7|, one can get solutions which have 

the minimum action within each twist sector. By taking pi = Ni, p' { = N 2 , 
s i,2 = s\ 2 = 0, and s\ N — S3 Ni = 1, one gets Q = 1/N. Other choices are 
easily extracted from Eq. 

8.4 General non-abelian minimum action solutions 

Concerning non-abelian solutions in arbitrary sizes and topological charges 
not much is known analytically. Mathematicians have given existence proofs 
of some of these configurations JJIJ. However, a good deal is known about 
the form and properties of these solutions by numerical methods. We will 
in what follows give a brief description of these results without entering too 
much into the technical details of the method employed. For that we refer 
the reader to the original papers pH], 



8.4.1 SU(2) 

The numerical work in this case has concentrated in studying the minimum 
action solutions for the most spatially-symmetric choice of twist in — k — 
(1, 1, 1) and size (7i = l 2 = h = l s ). In this case, n(n^ v ) = | and the action 
of all configurations must be larger than 47r 2 . The numerical method used 



by the authors of Ref. |3T[ is based on the lattice gauge theory formulation 
of Wilson [|l| . This is a gauge- invariant discretization of the Yang-Mills the- 
ory. As mentioned previously, putting the system on a box with periodic 
boundary conditions amounts to the discretization of a toroidal base space. 
Furthermore, it is possible to give in this case a lattice counterpart of the 



twisted boundary conditions [15|, |i9| . Although the lattice formulation de- 
stroys some of the topological features of the continuum theory, this is not 
the case for twist. Then, the point is to find by numerical methods the 
configuration which minimizes a discretized version of the continuum action- 
functional. The authors have used different discretized versions, but most 
of their results are based on Wilson's discretized form. By employing a so- 
called classically improved lattice action, one makes the discretization errors 
smaller, but unfortunately the numerical techniques become more time con- 
suming. In order to obtain the minimum lattice action configuration, the 
authors of Ref. |31"| , |32|j used a local minimization method known as cool- 
ing ]34j , |35| . Finally, they extracted discretized estimates of the relevant 



quantities and studied their convergence towards the continuum limit as the 
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Figure 1: We plot the value of the lattice action S divided by 4tt 2 as a 
function of the square lattice spacing a 2 , for the SU(2) configuration which 

— * 

minimizes the action with twist m — k — (1,1, 1). The solid line is a linear 
fit to the data. 
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lattice becomes finer. 

The lattices are hypercubic and with sizes which were always symmetric 
in space (N% x N t ), where N s and N t are the number of points in the spatial 
and temporal directions. If one fixes the continuum unit of length by setting 
l s — 1, then the lattice spacing is given by a = 

The simplest quantities to study are global ones, such as the total ac- 
tion and topological charge. Fig. 1 shows the value of the Wilson action 
for different values of a and N t /N s — > l t /l s m 2. The data approaches 4tt 2 
linearly in 1/a 2 . The slope depends on the lattice action used, but the ex- 
trapolation can give very precise results. From the data of Fig. 1 one gets 
an extrapolated continuum action of 4tt 2 with a precision of 1 part in 10 5 , 
which is consistent with the size of the errors from other sources. At the 
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same time, the topological charge tends to Q = | with equivalent precision. 
This is a strong indication that one is indeed approaching a self-dual clas- 
sical configuration. This is confirmed by studying self-duality directly (See 



Refs. [[31], [32], [36| for details). The other quantities studied also show that the 
lattice minimum configuration approaches a continuum one. The accuracy 
of the approximation depends on the sizes, the quantity studied, the choice 
of discretized action and the choice of discretized estimate of the different 
continuum quantities involved. In any case, for the values studied, the nu- 
merical differences were at most of a few percent. In what follows, we will 
omit all reference to the method and concentrate on the results, taken as 
valid for the continuum minimum action configurations. 

The next quantity to study is the energy profile of the classical solution, 
defined as: 

£{t)= [ d 3 xTr(E 2 (t,x) +B 2 (t,x)) . (96) 



For l t /l s > 1 this function has a maximum at a given time value, which will 
be called the time-center of the configuration, and chosen to correspond to 
t = 0. The main features shown by the energy-profile are the following: 

• The curve approaches a well defined one when l t /l s — > oo. 

• The curve is symmetrical about the time-center (t = 0). 

• For large values of \t\, £{t) goes to exponentially. 

• The curves for finite lt/l s > 1 differ from the infinite k/l s one mostly 
at the tails. Actually, the shape of the energy profile is well described 
by a periodization of the l t /l s — > oo curve. 

• The same curve is obtained within errors if one takes twice the electric 
part (or magnetic part) of £(t), in agreement with self-duality. 

• To show the quantitive aspect of the function, we show in Fig. 2, the 
shape of a function which fits the numerical data for the lt/l s oo 
curve, with a precision of the width of the line. The function has only 
4 parameters and its form is 

F(t) 



0.0113 + 0.00323 cosh(8M5t) + 0.0459t 5 
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Figure 2: We show the energy profile of the SU(2) minimal configuration 
with m — k — (1,1,1). 



Fig. 2 
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Notice that the coefficient of the exponential decay is roughly equal to 2%\/2 pa 
8.88, which is the minimum curvature of the classical potential. 

A qualitative conclusion is that the l t /l s — > oo solution can be consid- 
ered an instanton-like solution, which tunnels between 2 classical zero-energy 
states. The exponential decay at oo tells us that the system is not scale- 
invariant, being broken by the spatial size of the torus. We will refer to this 
object as Q = | instanton. 

In Refs. [[H], [32|, [36|] the authors investigate in some detail the form of 



several physical quantities for this configuration. The results are obtained 
for a ratio N t /N s = l t /l s pa 2 or higher, which except at the tails is already 
very close to oo. The field tensor F^x), the straightline Polyakov loops and 
the vector potential itself are studied. The latter is obtained in the A = ; 
Ai(t = — oo) = gauge. The main results obtained are: 

• The configuration action density has a maximum at a given point in 
space and time (center of the Q = \ instanton). The size of the object 
(Full width at half maximum) is a large fraction of the torus spatial 
size. 

• The configuration is invariant under spatial rotations belonging to the 
cubic group. This is the group which is left invariant by the boundary 
conditions. 

• The solution is highly non-abelian. At the vicinity of the Q = | in- 
stanton center, the solution behaves similarly to the ordinary BPST 
instanton. Hence, E\, E 2 and E 3 are orthogonal and of equal modu- 
lus in gauge space. In this sense it looks completely opposite to the 
constant field strength solutions (abelian and non-abelian), for which 
these quantities are parallel in Lie Algebra space. 

• Both the electric field(= magnetic field) and vector potential compo- 
nents can be nicely fitted with a few Fourier coefficients. 

• The space of solutions of this type {Q = \ for this twist) depends 
on 4 parameters, which can be taken as the coordinates of the center 
of the configuration. This, as will be commented later, is what the 
index theorem predicts for this case. In addition, there is a discrete 
degeneracy, about which the index theorem can say nothing. Indeed, 
there are actually 8 gauge inequivalent families of instantons. They 
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differ, for example, by the sign of the straight Polyakov loops passing 
through the Q = \ instanton center. 

8.4.2 SU(N) 

In Ref. B71 a similar study is performed for other SU(N) groups. Again, 
the study focuses in a spatially symmetric situation: h = h = h = h and 
fh = (1, 1, 1). Different groups are studied from N=3 to N=25. The reason 
for studying several groups will be more clear in the next section. Actually, 
there are two cases which are particularly interesting: SU (3), because it is the 
physically relevant case for QCD, and SU(oo), because there is a hope that 
the solution simplifies in that limit and hence a higher chance of achieving 
an analytical solution, which is one of the goals of the numerical study. The 
temporal twist is chosen in several ways, but such that \n{n^ u )\ = 1. Let 
us concentrate in this paragraph on the results for N equal or close to 3. 
Essentially, the conclusions are similar to the SU(2) case, except that now 
the cubic symmetry is replaced by the cyclic group x — > y —>■ z —>■ x. Another 
change is that at the Q = jt instanton center, the three spatial components 
of the electric field are no longer orthogonal. In this sense, as N grows the 
configuration looks more abelian at the center. The shape of the energy 
profile does not depend on the temporal twist (with \n(n^ u )\ = 1) for l t /l s 
large enough. The form of the energy profile does not change too much from 
one group to another, provided we scale the time variable by 1/N and profile 
itself by N 2 . The resulting curves N 2 £^(t/N) are shown in Fig. 3 for various 
groups. 

8.4.3 SU(oo) 

The large N limit is particularly interesting, since as 't Hooft 0] noticed, 
the configurations which have an action going like 1/N (unlike instantons), 
are not suppressed as N — - > oo in the path integral. Furthermore, as ex- 
plained previously there is a bigger hope to be able to obtain an analytical 
solution in this case. Unfortunately, the approach towards iV — > oo in some 
quantities is fairly slow. Exploring larger groups numerically is limited by 
computer memory and time. This slow tendency shows up, for example, in 
the behaviour of the energy profile. For N « 3, as mentioned in the previous 
paragraph, its behaviour suggests that there could be a well-defined limit of 
the function N 2 £ N {t/N) as iV — > oo. However, although numerically quite 
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Figure 3: Comparison of the energy profiles for different SU(N) groups. The 
profiles are scaled as shown in the figure N 2 £ N {t/N). 



Fig. 3 
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similar, there is no clear way to extrapolate to SU (oo). Actually for groups of 
the order of A^ = 19 the profile shows a double maximum structure, instead 
of the single maximum one. 

There are, nevertheless, several features which quite consistently seem to 
become more exact as A" — > oo. One of the most remarkable is that gauge 



invariant quantities tend to be constant in space [33]. The final conclusions 
of this analysis are currently under study (Ref.pTfl). 



8.5 Other results on classical solutions 

In this section, we would like to comment about the result obtained by 



P. van Baal and P. Braam |38| concerning the non-existence of Q = 1 in- 
stanton solutions on the torus without twist (See also [^9| ). A recent less 
technical discussion of this result is done in Ref . EDU . The result is based on 



Nahm's transformation [f|l], This transformation takes an U(N) self-dual 
gauge field A M (x) on the torus with topological charge Q, and gives rise to a 
new self-dual U(Q) gauge field living in the dual torus, and whose topologi- 
cal charge is N. If we start with an SU(N) gauge field, the transform is an 
SU(Q) gauge field. As a consequence, it is clear that there cannot exist an 
SU(N) self-dual gauge field with unit topological charge, since its Nahm's 
transform would make no sense. 

For completeness, let us sketch the definition of the Nahm's transform 
for the gauge potential A jU (z). First, one introduces a whole family of gauge 
inequivalent self-dual gauge fields as follows: 

A„{x, z) = A M (z) + 27r Zf J , (97) 

where < < are real numbers. They label the points of the dual torus 
of size -^x . . .x j^. This is so, because there exist a gauge transformation that 
maps A^(x, z v ) into A M (x, z v + l/l u ). The next step is to consider the Dirac 
operators corresponding to this family of gauge potentials. It is crucial that 
these operators have exactly Q positive chirallity (ip^(x, z) for i=l . . . Q) and 
no negative chirallity solutions (the difference is fixed by the index theorem). 
Then the Nahm's transform is given by: 

dxiP Hr) (x,z)—^\x,z). (98) 
The proof that A l Hz) is self-dual, follows very closely the ADHM construc- 



tion 43, 44 . 
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It is tempting to think that the Nahm's transform could allow the con- 
struction of new self-dual solutions, starting from known ones. However, it 
is known that the Nahm's transform of a constant field strength solution 
is another constant field strength solution. An explicit example is given in 
Ref.pl. 



9 Fluctuations around classical solutions 

In this section we will very briefly comment a few things about fluctuations 
around classical solutions. The first point we want to mention concerns zero 
modes. If we have a classical solution and we consider a small perturbation 
around it, the variation of the action is to first order equal to zero. To second 
order, we have a quadratic form given by an operator whose spectrum is the 
spectrum of fluctuations. For a minimum action solution all eigenvalues are 
positive or zero. The corresponding zero-eigenvalue eigenvectors are called 
zero-modes. If the classical solution depends on several parameters, the vari- 
ation with respect to each parameter gives rise to a zero-mode. Hence, there 
are at least as many independent zero modes as parameters of which the 
solution depends. For a gauge theory, there are always zero-modes present, 
associated with gauge transformations of the solution. Hence, what one is 
really interested in, is in those zero-modes, giving rise to parameters, which 
are not associated to gauge transformations. The space of these parameters 
is what is called the moduli space. The number of independent zero-modes 
not associated with gauge transformations gives its dimensionality. 

For the case of self-dual solutions, it was found some time ago H5I, [46] , that 
the Atiyah-Singer index theorem can tell us the dimensionality of the moduli 
space d, even if we do not know the analytic expression of the solutions. 
For ordinary instantons this is well-known: SU (2) self-dual solutions on the 
sphere are known to depend on d = 8\Q \ — 3 parameters (Q is the topological 
charge). This was known before the ADHM method gave us a recipe to 
construct them. However, the result depends both on the gauge group and 
on the base-manifold (See Refs j|7|, f48|). In particular, for the torus T 4 the 
result is d = 4ca(G) \Q\, where ca(G) is the quadratic Casimir in the adjoint 
representation of the group G. Hence, SU(N) self-dual solutions depend on 
d = 4 N \ Q\ parameters. 

As a consequence, of the previous formula, the Q = ■h instantons which 
appear in twisted SU(N) theories, depend on 4 parameters. These param- 
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eters are associated to the translational zero-modes in 4 dimensional space. 
There is no zero-mode associated with scale invariance, as for BPST instan- 
tons, since this symmetry is broken by the torus size. Actually, the form of 
d is consistent with the interpretation of the space of self-dual solutions, as 
a liquid of Q — instantons done in Ref. (T(| . 



An analysis of the fluctuation spectrum around the constant field-strength 
solutions of subsection 7.2 has been carried by P. van Baal He concluded 
that the only solutions which are stable are the self-dual (or anti-self-dual) 
ones. The eigenfunctions are given in terms of generalized Riemann theta 
functions. 



10 Concluding remarks 

In the previous pages, we have reviewed some of the most important results 
obtained for classical Yang-Mills fields on the torus, many of which turn out 
to be relevant when one studies the quantum mechanical system. The second 
part of my lectures at Peniscola, covered this aspect. However, the written- 
up version of this part is not yet ready and will appear later. Unfortunately, 
a good deal of our motivation for plunging into the study of the classical 
system lies precisely in the second part. We hope this will not discourage 
the interested reader. 

To conclude, I will just mention that some topics have been left out due 
to lack of space and time. For example, we have not included the results 



obtained about sphalerons in recent years [p0|, |51fl . In any case, I hope that 
the present review will turn to be useful for those interested in studying the 
subject. 
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